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ABSTRACT 

We predict the redshift of the first observable (i.e., in our past Hght cone) Gamma 
Ray Burst (GRB) and calculate the GRB-rate redshift distribution of the Population 
III stars at very early times (z = 20 — 60). Using the last 2 years of data from Swift 
we place an upper limit on the efficiency {tjcrb) of GRB production per solar mass 
from the first generation of stars. We find that the first observable GRB is most likely 
to have formed at redshift 60. The observed rate of extremely high redshift GRBs 
(XRGs) is a subset of a group of 15 long GRBs per year, with no associated redshift 
and no optical afterglow counterparts, detected by Swift. Taking this maximal rate we 
get that riGRB < 1-1 10~^ GRBs per solar mass in stars. A more realistic evaluation, 
e.g., taking a subgroup of 5% of the total sample of Swift gives an upper limit of 
VGRB < 3.2 10-5 GRBs per solar mass. 

Key vifords: galaxies :high-redshift - cosmologyitheory ~ GRB 



1 INTRODUCTION 

Gamma Ray Bursts (GRBs) are the bri ghtest known events 
■ since the Big Bang (see reviews by: IPiranI |2005| . l2000l : 
, iMeszaroi |2002| ) . Therefore they offer a wonderful prospect 
' to explore the evolution of the Universe ever since stars be- 
gan to fo rm. The f irst generatio n of stars is ex pected to be 
. massive (|Abel et al . 2002; Bro mm. Coppi. fc L arson 2002), 
' and at least som e of them should f orm a GRB as they end 
their lives (e.g., iHeger et al.ll200^ ). Since the luminosities 
of galaxies and quasars decline with redshift, high redshift 
GRBs are expected to be easier to observe than their host 
galaxies. Consequently, they can offer a signific ant probe of 
the early universe fe.g.. iLamb fc Reichartll2000l ). For exam- 
ple, since the GRB afterglows fade only slowly with red- 
shift they offer a source for stu dying the cosmic reionization 
(IBarkana fc Loebll2004l2006l ). In particular, GRBs forma- 
tion history is expected to fol l ow the star formation ev olu- 
tion l|Blain fc Nataraianllioool : iPorciani fc Madaull200ll ). 

The first stars in the Universe namely, population 
III (POP III), formed at very high redshift (z < 66, 
iNaoz. Noter fc Barkanal bood ). Since metals are absent in 
the pre-stellar universe, the earliest available coolant is 
molecular hydrogen (H2). Thus the minimum halo mass that 
can form a star is essentially set by requiring the infalling 
gas to reach a temperature > 1000 K required for exciting 
the rotational a nd vibrational states of m olecular hydrogen 
llTegmarklll997l). Numerical simulations ( Abel et al.1 l2002l: 
iFuUer fc Couchman .1120001 : lYoshida. et al.|[2003l : iReed et aH 



I2005I : n romm. Coppi. fc LarsonI I2OO2D give a more accu- 
rate constraint and require a minim um circular velocity 
K ~ 4.5 km/s, where K = y/GM/R in terms of the halo 
virial radius R. These simulations include gravity, hydro- 
dynamic and chemical processes in the primordial gas, and 
show that the first star formed within a galactic halo of 
~ 10* Af© in total mass. The predicted mass of this star is 
quite heavy and should exceed IOOMq (jYoshida et al.ll2006l : 
iGao et al.ll2007l ). 

The radiation from these first stars is expected to even- 
tually dissociate all the H2 in the intergalactic medium, lead- 
ing to the domination of a second generation of larger galax- 
ies where the gas coo ls via radiative transitions in atomic 
hydrogen and helium (|Haiman et al.|[l997 ). Atomic cooling 
occurs in halos with Vc > 16.5 km/s, in which the infalling 
gas is heated above 10,000 K and is ionized. 

iBromm fc Lo^ { 20061 ) calculated the star formation 
evolution for population III while assuming only atomic 
cooling in order to derive the GRB rate. They assumed 
that the efficiency of GRB production (tigrb) is ~ 10"^ 
GRBs per solar mass in stars, and found that about 
10% of all bursts detected by Swift should be generated 
from redshift ^ 5. Unlike the se writers and GRB rate 
analys is done by others (e . g .. iGuetta. Piran. fc WaxmanI 
'2003 : iGuetta fc PiranI |2005|, l2007l: i Bromm fc Loebl l200a 
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therein), we concentrate on a much higher redshift regime 
{z = 20 — 60) and we estimate the upper bound of the effi- 
ciency of extremely high redshift GRBs (XRGs) production. 
The article is structured as follow: we begin by presenting 
a simple star formation evolution (Section [2]| , relevant for 
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these redshifts. Our analysis is presented in Section [3] In 
Section [3.11 we calculate the redshift of the first observable 
GRB. Section [3.21 is dedicated to the calculation of the red- 
shift distribution of the XRG rate, in order to place an upper 
limit on rjcRB ■ We conclude with a discussion of our results 
(Section II]). 

Our calculations are made in a ACDM universe, in- 
cluding dark matter, baryons, radiation, and a cosmological 
constant. We assume cosmological parameters matching the 
three ye ar WMAP data tog ether with weak lensing obser- 
vations (|Spergel et aLlbOod ). i.e., erg = 0.826, = 0.299, 
= 0.74, and Qi = 0.0478. 



2 ESTIMATION OF THE HIGH REDSHIFT 
STAR FORMATION RATE 

We adopt a simple model for the star formation history at 
high redshift (z ~ 70 — 20). This modle assumes that at 
these early times most of the stars formed out of newly ac- 
creted gas during mergers. We follow the linear and no n lin- 
ear perturbation growth from iNaoz fc Barkanal (|2005l ') and 
iNaoz. Noter fc Barkanal l|2006l ) in order to derive the frac- 
tion ofjnassjrijialos. 

iNaoz fc Barkanal l|2005h showed that the baryon sound 
speed varies spatially, so that the baryon temperature and 
linear density fluctuations must be tracked separately. In 
the non-linear regime, it is necessary to calculate correctly 
the linearly extrap olated overdensity Sc wh i ch m arks the 
time of the collapse. iNaoz. Noter fc Barkanal l|2006l ) showed 
that the value of Sc is lower at high redshift than the 
classical value (1.686) and varies with time (see for refer- 
ence fig. 6 in iNaoz fc Barkanal l2007h . Defining the mass 
varian ce 5* = a^{M,z), we calculate the ISheth fc TormenI 
l|200ll ) mass function, which fits simulations and includes 
non-spherical effects on the collapse. The function fsT is 
the fraction of mass associated with halos of mass M: 



fsT{5c,S)=A' 



sy 2tt 



1 + 



exp 



(1) 



where u = Sc/y/S. We use best-fi t parameters a' — 0.75 and 
q' = 0.3 (jSheth fc TormenI I2OO2I ). and ensure normalization 
to unity by taking A' = 0.322. We apply this formula with 
5c{z) and a^{M,z) as the arguments. 

As mentioned, the only coolant available for these high 
redshift halos is cooling via H2 cooling. Thus, we find the 
fraction of mass in halos with mass larger than the minimum 
H2 cooling mass to be 



F{> Mrmn,H,{z)) = / /sT (5c (z) , S)dS 



(2) 



Thus the star formation rate (SFR) (based on H2 cooling) 
is simply: 



SFR - 



dF Q.b 



(3) 



where po is the com oving matter den s ity, a nd e is the star 
formation efficiency. iBarkana fc Loebl l|2000t ) calculated the 
star formation history for lower redshifts using a more com- 
pHcated calcula tion that included addi tional merger-induced 
star formation. iBromm fc Loebl (|2006l ) have calculated sep- 
arately the star formation rate for the different star popula- 
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Figure 1. Redshift evolution of SFR in halos larger than the 
minimum H2 cooling mass, in Mq yr~^ per comoving Mpc'^. 
This is the prediction from the model presented in equation |3] 



tions, i.e., POP I/II and POP III. In both papers the star for- 
mation efficiency was chosen to be 10% independent of red- 
shift. This yields a rough agreement betwe en the SFR at low 
redshift and the observations (see fig. 1 in lBarkana fc Loebl 
[20001 . Following these previous analyses we set e = 10%. 

In figure [1] we plot the SFR as a function of 1 + 2 as 
derived from eq. |3| We note that th is resu l t is v ery close 
to that obtained by iBromm fc Loebl (|2002l , |2006| ) for the 
star formation of POP III. However our result predicts a 
higher rate because we consider star formation at higher 
redshifts. In practice our result is between their most opti- 
mistic (high efficiency) estimate and their most pessimistic 
case. In addition our simple model does not include the ra- 
diation feedback which lowers the contribution of molecular 
hydr ogen cooling compared to atomic cooling (s ee for refer- 
ence |Haimaner2]|l993; llarfena^L^ Therefore 
we restrict ourselves to al l redshifts ^ 2 (note that the red- 
shift range considered bv lBromm fc Loe b (2006) is z ^ 30). 
Since we do not account for reionization effects at all we 
compare our r esults to that of a lat e reionization in (figure 
1 top panel in iBromm fc Loeblbood ). 



3 CALCULATIONS AND RESULTS 

In order to produce a GRB the progenitor must ful- 
fill 3 requirements (see Zhang. Wooslc v. fc Heged |2004| : 
iPetrovic et al ] |2005I :Ib romm fc Loehi2006i . for more detailed 
explanation). First, the star must be massive enough to re- 
sult in a black hole. Second, the star must lose its hydrogen 
envelope in order to allow the relativistic jet to pen etrate 
through the surface (|Zhang. Wooslev. fc Hegen 120041 ). and 
finally an accretion disk must be able to form. For that the 
star must maintain enough angular momentum. When eval- 
uating the efficiency of producing GRBs from POP III stars 
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there are many uncertainties, such as the properties of the 
progenitor. iBromm fc Loebl l|200d ). for example, claim that 
in order to produce a GRB, the progenitor must have a close 
binary companion that helps in stripping the hydrogen enve- 
lope and retaining suf ficient angular momentum in th e col- 
lapsing core. However. iFrver. Wooslev. fc Heeerl ()200ll ) show 
that massive POP III stars can produce pair-instability su- 
pernovae that can strip off their hydrog en envelope and re - 
sult in a black hole accretion disk system. [Heger et al.l(l2003l ) 
deploy a variety of masses that can produce GRBs from sin- 
gle POP III stars. Whether or not binary stars are a strict 
requirement for GRB production is still debatable. However, 
our analysis is independent of that question. We only assume 
that GRBs in this high redshift regime do exist, and use ob- 
servational constrains to extract the efficiency of producing 
them. 

The most naive assumption is that each star (or binary 
system) at the relevant high redshift range produces a GRB 
and thus tigrb ~ 1/100 GRBs per solar mass (or 1/200 
assuming binary progenitors) . We derive further restrictions 
on riGRB by normalizing the GRB rate with the potential 
XRG rate observed by Swift (see text below for details). 



3.1 The First Observable GRB 

iNaoz. Noter fc Barkanal (|2006f ) have found that the first star 
most likely formed ai z = 65.8. As mentioned above, this 
star is expected to have a mass > 100 Mg l|Abel et al.l 
|2002| ). thus it (as well as the stars that follow it) is a suit- 
able candidate to produce a GRB (|Heger et al.ll2003l ). We 
assume that GRBs at these redshifts are bright enough to 
be observed today (estimates of the limiting luminosity are 
given below using a more compl ex model) , and ad opt a fidu- 
cial jet opening angle of 5° (e.g.. lFrail et al.ll200ll ). This lim- 
ited opening angle means that we see only two out of 500 
GRBs. The probabilit y to observe the first GRB on t he sky 
is computed following iNaoz. Noter fc Barkanal (|2006l ) anal- 
ysis for the formation of the first star. Notice that since 
the lifetime of the first stars is insignificant compared to the 
Hubble time, we can neglect the delay time between the star 
formation and explosion. We find that most likely the first 
observable GRB is formed a.t z — 59.9, with a l-cr (68%) 
range of 2 = 58.9-61.4 and a 2-a (95%) range oi z = 58.1- 
63.5. Note that the most probable redshift of the first GRB 
is only weakly sensitive to changing the opening angle. For 
example an isotropic GRB will be seen at the redshift of the 
first observable star. 20° we find that Considering an open- 
ing angle of 5° we place an upper redshift limit for observing 
GRBs. Combining with the limitation of our SFR model we 
concentrate on the redshift range of: 20 — 60. 

3.2 GRB production efficiency and rate 

From equation ([3]) we can evaluate the observable GRB rate: 

Rgrb{z) = r)GRB f^^ TT- / (t>{L)d\ogL , (4) 

1 + z az 

where V is the comoving volume and the factor (1 -I- z)~^ 
accounts for the cosmological time dilation, tigrb was de- 
fined earlier as the GRB efficiency in units of GRBs per solar 



mass, and 4>{L) is the GRB luminosity function. The lim- 
iting luminosity at a given redshift, Lmin, depends on the 
flux sensitivity threshold of the detector [fum), and can be 
written as: 

T ,U \ C{El{l + z),E2{l + z)) Y^ 

Lmin ^ 4TTdL[z)tli„,< ^^^^ > , (5) 

where dL{z) is the luminosity distance and C(E1,E2) is the 
total integrated luminosity for observation with Swift be- 
tween _E1 = 50 and E2 = 300 ke\{^ (for specific expla nation 
see for reference [Cuetta fc Piranll2007l : [Schmidtlll999^ . 

Assuming an average photon number distribution of a 
single power la w, eq. ((SJ can be written as (see for reference: 
ISchmidt|[200ll ): 

Lm^n = 4Tvdl{z)fum{l + z)"'^ + "' . (6) 

where a is the average photon spectral index over the en- 
ergy range (i.e., dN/dE oc _E", where A'^ is the number of 
photons). A typical GRB photon number distribution can 
be describ ed as a broken po wer law, also known as a Band 
spectrum ([Band et al.lll993l ). The low energy regime has a 
spectral slope of a ~ — 0.8, while the high energy regime 
has a slope of /3 ~ —2.5 jPreece et al.|[200ol ). The peak en- 
ergy {Ep) falls in the range of 100 keV to a few MeV, which 
renders it observable by Swift. 

Several groups (e.g., lYonetoku et al] 12004 : 
iGhirlanda et all l2005l : iLamb. Donaghv. fc Grazianil 12005 ) 
using various samples taken from BATSE, BeppoSax, 
Integral and HETE2, found that Ej, ~ Llf^, where Lisa is 
the GRB isotropic equivalent luminosity (either average 
or peak luminosity). Each burst that has an Ep which 
falls within our observation window, can be fitted with an 
average photon spectral slope, —0.8 ^ a ^ — 2.5. Because 
Ep is the break energy in the power law and bursts with 
higher luminosity have higher Ep, the average observed 
slope of such bu rsts should be higher (i.e., shallower) as well. 
ISchmidtl l|200ll ) found for a sample of simulated long GRBs 
a median value of a = —1.6. In addition he found a positive 
correlation between a and the GRB luminosity, as expected 
from the Ep-Liso relation. For example a typical XRG, at 
z = 20, with an intrinsic luminosity of 10^^ erg sec~^ has 
an Bp ~ 1 MeV which will be redshifted to ~ 50 keV. The 
photon spectrum within the BAT observational window 
will have an average spectral slope of a ~ —1.6. 

We use this single power law to get a simplified k- 
correction for our bursts, where we expect most XRG spec- 
tral slopes to fall between —1.6 and —1, since their luminos- 
ity is in the high tail of the luminosity function. We adopt 
a limiting case of a = — 1, according to the Schmidt (200l|) 
hardness luminosity correlatior(f]. We also address the differ- 
ence that adopting the median value (—1.6) would make in 
our results. 

iLamb fc ReichartI (|2000l ) assumed that the Swift thresh- 
old limit would be fum = 0.04 photons cm~^ sec~^, and 



^ The actual observed energy band of Swift is : 15 — 150 keV, while 
our calculation follows FCuetta fc PiranI | |2007|) . who extrapolated 
to the BATSE energy band. 

^ This is of course the more interesting case since it favors higher 
luminosity bursts and suits the redshift regime in this paper. 
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that the spectral photon index is —1. For an observed pho- 
ton number flux with an observed redshift they found the 
luminosity of the burst as well as the maximum observable 
redshift. For example, for burst GRB971214 they found a 
peak luminosity of 6.4 x 10^® sec~^, and showed that a burst 
with this luminosity can be seen up to redshift ~ 70. Re- 
peating this exercise with the 2-year Swift threshold limit 
fiim = 0.2 photons cm~'^ sec~^ we find that a burst with 
this peak luminosity can be seen at least up to z ~ 21, while 
adopting a — —1.6 yields a maximum observable redshift of 
7.8. Thus, the maximum observed redshift for a given lumi- 
nosity is dependent on the spectral index. 

In figure [2] we plot the maximum observable redshift of 
a burst as a function of the spectral index, for three different 
representative luminosities, from bottom to top: 10""^, 10^^ 
and 10^^ ergs sec~^. The figure is plotted for fum = 0.2 
photons cm~'^ sec~^, and it can be seen that the maximum 
observable redshift of a fixed luminosity GRB is strongly 
dependent on the spectral index. Note that high luminosity 
bursts are observable to higher redshift, especially since they 
tend to have high a. In the figure we indicate the limiting 
value and the median value of a , —1 and —1.6, respectively. 
Thus, high luminosity bursts (~ 10^^ ergs sec~^) with a ~ 
— 1 can be expected to be observed up to 2; ~ 80, while 
for the median value of the spectral index —1.6 the bursts 
are expected to be observed up to z; ~ 18. Thus observing 
GRBs from a higher redshift (while assuming a = —1.6) 
suggests much more luminous bursts than ~ 10^^ ergs sec~^. 
Of course this can also mean (when placing an upper limit 
on the luminosity) that there are no observable extremely 
high redshift GRBs with the current Swift detection limit. 
Therefore for Swift to detect XRGs, their average spectral 
slope should be closer t o —1, as shown in fig. [H 

We assume th e iGuetta. Piran. fc WaxmanI (|2005l ) 
and Guetta fc PiranI (|2007l ') luminosity function based on 
ISchmidtl l|200ll ). This luminosity function has a broken 
power law with cutoff luminosities of L*/Ai and A2L* . 



<f>{L) 



Co 



Ai 

L* < L < A2L* 



(7) 



where (j>{L) is normalized so that the integral over the lumi- 
nosity function e quals unity through co. We have used the 
iGuetta fc PiranI (|2007l ') best fitted values with Ai , A2 = 100, 
L* = 6.5 X lO'^^ erg sec"\ a = -0.2 and b = -1.7, marked 
in their paper as model (v). This model favors high lumi- 
nosity bursts, which are the interesting case here. We also 
address their (vi) model with b — —2 which has fewer high 
luminosity bursts (see text below). 

We can now return to eq. Q and calculate the XRG 
rate from POP III stars at z = 60 — 20, and use that to 
place an upper limit of the efficiency of GRB production 
from these stars. We ask what is the efficiency of GRB pro- 
duction assuming that Swift observes up to ~ 15 long GRBs 
per year originating from these high redshift POP III stars. 
This subsample includes GRBs with no associated redshift 
(there are about 60 per year out of 90 per year). In addition 
we exclude GRBs with Tgo shorter than 50 sec0, where Tgo 



A burst with T90 = 50 at 2: = 20 has a proper time duration of 
~ 2.5 sec which places it on the border between long and short 
GRBs. 
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Figure 2. The maximum detectable redshift by Swift as a func- 
tion of the spectral photon index a, using the Swift detection 
threshold fu^ = 0.2 photons cm~^ sec~^. We consider three 
representative luminosities. From bottom to top: 10^^,10^^ and 
10^'^ ergs sec~^. Als o shown are the median spectral index value 



from lSchmid^ bOOih 



-1.6 and the limiting case a ■■ 



defines the duration of the burst, since short GRBs most 
likely origin ate from much lower redshift (see for reference: 
lNakaill2006h . Finally we exclude GRBs with no optical coun- 
terparts for the following reason: consider a burst at 2 ~ 20, 
with an afterglow counterpart. The photons of the burst red- 
shift with the expansion of the Universe as they propagate. 
Wavelengths shorter than the Lyman limit will be absorbed 
by photoionizing hydrogen or helium atoms. Assuming that 
the Universe was completely reionized at 2 ~ 10 or later each 
2: = optical photon was at a wavelength shorter than the 
Lyman limit and thu s was absorbed in the p re-reionization 
era (see for reference: iBarkana fc Loebll200ll ). Therefore, no 
optical afterglow should be observed from XRGs, and we are 
left with about 15 bursts per year from the Swift sample. Of 
course assuming that all of this sample originates from high 
redshift is mostlikely an overestimate. Nevertheless, it gives 
an upper limit to the efficiency. Below we give more realistic 
estimates of the XRG rate detected by Swift. 

In figure [3] top panel, we depicted the predicted GRB 
rate redshift distribution for various efficiency normaliza- 
tions. We consider both cases of q = —1 (solid curves) 
as well as a = —1.6 (dashed curves). For the first case 
we consider XRGs that contribute (from top to bottom) 
15, 5 and 1 GRBs per year to the total sample detected by 
Swift. The GRB efficiencies from these values are: rjcRB = 
1.1 10"'',3.5 10"^ and 7.1 10"^ GRBs per solar mass, re- 
spectively. For the latter case {a = —1.6) we consider the 
same contributions, (top to bottom) 15, 5 and 1 GRBs per 
year from the relevant redshift. The resulting upper limit 
on the efficiency is r)GRB = 2.6 10"^, 8.7 10"" and 1.7 10"" 
GRBs per solar mass, respectively. Not presented in the fig- 
ure is the lower limit case of a — —2.2. We find that for 
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this spectral index Swift will be able to detect only a few 
GRBs per year. Comparing with the results above, in this 
case the efficiency needed to contribute one GRBs per year 
detectable by Swift is 4.2 10""^ GRBs per solar mass. 

As can be seen from the figure in order to generate the 
same contribution of GRBs per year, the efficiency must be 
different for different spectral indices. Adopting a = —1 gen- 
erates a lower efficiency for the same percentage of GRB con- 
tribution. This can be seen more clearly in the bottom panel 
of this figure, where we depict the GRB efficiency per solar 
mass (jigrb) for various fractional contribution of GRBs 
from redshifts 20 — 60 (out of ~ 90 GRBs per year detected 
by Swift). We consider the spectral index of —1 (solid curve) 
and also spectral index of —1.6 (dashed curve). For example, 
for a 1% contribution of GRBs detected by Swift per year 
from this redshift range, i.e., about one GRB per year, we 
find that a spectral index of —1 predicts about 2 orders of 
magnitude lower efficiency of GRB production than assum- 
ing a spectral index of —1.6. As mentioned before, the latter 
spectral index suggests much more luminous GRBs which 
is less reasonable. We also consider (in this^ bottom panel 
of figu re [3l) the parameters of model (vi) in iGuctta & P iranI 
l|2007l ) (doted curve), which is less favorable for high lu- 
minosity bursts (fe = —2). As can be expected this model 
results in a higher GRB efficiency since the luminosity func- 
tion doe s not favor h igh luminosity bursts. 

Br omm fc Loebl J2OO6,) found that about 10% of the 
GRBs detected by Swift each year should result from high 
redshift (^ 5) stars. However they assumed a lognormal 
luminosity function, and efficiency of ~ 10~^ GRBs per 
solar mass. Thus if this 10% claim is true, it is reason- 
able to assume that even less originate from a higher red- 
shift (^ 20) range. Suppose we assume a more reason- 
able upper limit for example, about 5% to the total sam- 
ple (about half of th eir range) . Using th eir luminosity 
function (see eq. 8 in iBromm fc Loebl l200 2h we find that 
for a = — 1, this 5% contribution (about five GRBs per 
year) results in an upper limit on the efficiency of: 1.8 10~® 
GRBs per solar mass, while assuming a = —1.6 gives 
an efficiency of 1.0 10~^ GRB s per solar mass. Adoptin; 
the Guetta. Piran. fc WaxmanI (|2005|) and iGuetta fc PTra . 
(I2OO7I ) luminositv function, (eg. ((7}) we have (as depicted in 
fig- 13) noRB = 3.2 10"^ and 7.8 lO""* GRBs per solar mass, 
for a — —1 and —1.6, respectively. 



4 DISCUSSION 

We have predicted that the redshift of the first observable 
GRB is most likely z ~ 60. Th is prediction is based on the 
iNaoz. Noter fc Barkanal (|2006f ) prediction of the redshift of 
the first observable star. 

We have calculated the GRB rate from extremely high 
redshift (z = 20 - 60) POP 111 stars. We used a simple SFR 
(based on H2 cooling) that includes the correct evolution of 
fluctuations. We assume that at these early times most of 
the stars formed out of newly accreted gas during mergers. 
We find that our SFR is consistent wi th previous analysis 
at lower redshift (|Bromm fc Loebll2006l ). However our result 
predicts a higher rate (see fig. [TJ because we consider star 
formation at higher redshifts. 

Using this SFR we evaluated the GRB rate redshift dis- 
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Figure 3. Top panel: Redshift distribution of GRBs for various 
POP III contribution to the total Swift detection sample. We con- 
sider a spectral index of —1 (solid curves) and also a spectral index 
of —1.6 (dashed curves). For the first case {a = —1) we consider 
a contribution of (top to bottom) 15, 5 and 1 GRBs per year from 
redshifts of 20 - 60, which results in tjgrb = 1-1 10-*,3.5 lO''^ 
and 7.1 10~® GRBs per solar mass respectively. For the second 
case {a = —1.6) we consider (top to bottom) 15,5 and 1 GRBs 
per year from the relevant redshifts. The resulting upper limit on 
the efficiency is r]GRB = 2.6 10-^,8.7 lO""* and 1.7 lO""* GRBs 
per solar mass, respectively. Lower panel: tigrb for various frac- 
tional contribution of GRBs out of ~ 90 GRBs per year detected 
by Swift, from redshifts 20 — 60. We consider the spectral index of 
— 1 (solid curve) and also spectral index of —1.6 (dashed curve). 
F or this latter case we als o consider the parameters of model (vi) 
in lGuetta & PiranI ||2007| ) (dotted curve), which is less favorable 
for high luminosity bursts {b = —2). 



tribution (eq. (jU). We adopted [Guetta fc PiranI (|2005l l lu- 
minosity function, and performed a simple k-correction. We 
showed that the maximum redshift to which GRBs can be 
observed depends on their luminosity and on the spectral 
index (see fig. [2ll. Using a l so the Ep-Liso relation and, as 
supported by the ISchmidtl (|200ll ) luminosity hardness cor- 
relation, we adopted a = —1. This value enables high lu- 
minosity bursts to be seen to larger distances, i.e., higher 
redshift. 

Using the Swift detection rate we have placed an up- 
per limit on the efficiency of GRB production from POP 111 
stars. For various fractional contribution to the total Swift 
sample we find different efficiencies (fig. [Sjl . The safest up- 
per limit we set is: 1.1 10^* GRBs per solar mass. This is 
a result from a maximal contribution of XRGs to the Swift 
sample. This sub-sample has 15 very long (Tgo > 50 sec) 
GRBs per year with no associated redshift and no optical 
afterglow counterparts (about ~ 17%, of the total Swift de- 
tection sample). We also find a more realistic estimate, that 
the upper limit of GRB efficiency from POP III stars is 
3.2 10~^ GRBs per solar mass. This efficiency results in a 
high redshift GRB contribution of 5% of the Swift sample. 
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Taking the median value a = —1.6 we find that the upper 
limit on tigrb is 7.8 10~^ GRBs per solar mass. 

We caution that we have extrapolated the properties of 
the observed GRBs to much higher redshift. An additional, 
important caveat to our model is the strong dependence on 
the spectral index a. We point out that given the Ep-Liso re- 
lation it is reasonable to assume that a ^ — 1.6. The E^-Liso 



relation is closly relat ed to the Amati et al. (2002) relation. 
Sever al authors (e.g.. iBand fc Preecai2005. : Nakar fc PiranI 
I2OO5I ) claimed that this relation may be a boundary curve 
created by the brightest and softest bursts, as a result of ob- 
servational biases. They suggested that many bu rsts should 
be dimmer and harder (though see lAma±ill2006l ). This sup- 
ports our argument that XRGs should have high a. For com- 
pleteness, though, we have shown that if a = —2.2, Swift can 
detect only a few GRBs per year from redshift > 20, and 
for one GRB per year we find a high upper limit on the 
efficiency of 4.2 10"'^ GRBs per solar mass. 
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